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Abstract 

We study the transport properties of model networks such as scale-free and Erdos- 
Renyi networks as well as a real network. We consider the conductance G between 
two arbitrarily chosen nodes where each link has the same unit resistance. Our 
theoretical analysis for scale-free networks predicts a broad range of values of G, 
with a power-law tail distribution <I>sf(G) ~ G"^'^, where go = 2A — 1, and A is 
the decay exponent for the scale-free network degree distribution. We confirm our 
predictions by large scale simulations. The power-law tail in ^sfCG*) leads to large 
values of G, thereby significantly improving the transport in scale-free networks, 
compared to Erdos-Renyi networks where the tail of the conductivity distribution 
decays exponentially. We develop a simple physical picture of the transport to ac- 
count for the results. We study another model for transport, the max-flow model, 
where conductance is defined as the number of link-independent paths between the 
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two nodes, and find that a similar picture holds. The effects of distance on the 
value of conductance are considered for both models, and some differences emerge. 
We then extend our study to the case of multiple sources, where the transport is 
define between two groups of nodes. We find a fundamental difference between the 
two forms of flow when considering the quality of the transport with respect to the 
number of sources, and find an optimal number of sources, or users, for the max-fiow 
case. A qualitative (and partially quantitative) explanation is also given. 
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1 Introduction 

Transport in many random structures is "anomalous," i.e., fundamentally dif- 
ferent from that in regular space [1,2,3]. The anomaly is due to the random 
substrate on which transport is constrained to take place. Random structures 
are found in many places in the real world, from oil reservoirs to the Internet, 
making the study of anomalous transport properties a far-reaching field. In 
this problem, it is paramount to relate the structural properties of the medium 
with the transport properties. 

An important and recent example of random substrates is that of complex 
networks. Research on this topic has uncovered their importance for real- 
world problems as diverse as the World Wide Web and the Internet to cellular 
networks and sexual-partner networks [4]. Networks describe also economic 
systems such as financial markets and banks systems [5,6,7]. Transport of 
goods and information in such networks is of much interest. 
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Two distinct models describe the two limiting cases for the structure of the 
complex networks. The first of these is the classic Erdos-Renyi model of ran- 
dom networks [8], for which sites are connected with a link with probability 
p and are disconnected (no link) with probability 1 — p (see Fig. 1(a)). In 
this case the degree distribution P(/c), the probability of a node to have k 
connections, is a Poisson 



where k = J2'kLikP{k) is the average degree of the network. Mathematicians 
discovered critical phenomena through this model. For instance, just as in 
percolation on lattices, there is a critical value p = Pc above which the largest 
connected component of the network has a mass that scales with the system 
size N, but below Pc, there are only small clusters of the order of log N. At 
p = Pc, the size of the largest cluster is of order of iV^/^. Another characteristic 
of an Erdos-Renyi network is its "small-world" property which means that 
the average distance d (or diameter) between all pairs of nodes of the network 
scales as logiV [9]. The other model, recently identified as the characterizing 
topological structure of many real world systems, is the Barabasi- Albert scale- 
free network and its extensions [10,11,12], characterized by a scale-free degree 
distribution: 



The cutoff value /cmin represents the minimum allowed value of k on the net- 
work (/Cmin — 2 here), and k^nax = ^minA^^^^'^~^\ the typical maximum degree 
of a network with N nodes [13,14]. The scale-free feature allows a network 

to have some nodes with a large number of links ("hubs"), unlike the case 




(1) 



P{k) ~ k 



-A 




(2) 



3 



for the Erdos-Renyi model of random networks [8,9] (See Fig. 1(b)). Scale- 
free networks with A > 3 have d ~ logA^, while for 2 < A < 3 they are 
"ultra-small- world" since the diameter scales as d ~ log log [13] . 

Here we extend our recent study of transport in complex networks [15,16]. We 
find that for scale-free networks with A > 2, transport properties characterized 
by conductance display a power-law tail distribution that is related to the 
degree distribution P{k). The origin of this power-law tail is due to pairs of 
nodes of high degree which have high conductance. Thus, transport in scale- 
free networks is better because of the presence of large degree nodes (hubs) 
that carry much of the traffic, whereas Erdos-Renyi networks lack hubs and 
the transport properties are controlled mainly by the average degree k [9,17]. 
We present a simple physical picture of the transport properties and test it 
through simulations. We also study a form of frictionless transport, in which 
transport is measured by the number of independent paths between source 
and destination. These later results are in part similar to those in [18]. We 
test our findings on a real network, a recent map of the Internet. In addition, 
we study the properties of the transport where several sources and sinks are 
involved. We find a principal difference between the two forms of transport 
mentioned above, and find an optimal number of sources in the frictionless 
case. 

The paper is structured as follows. Section 2 concentrates on the numerical 
calculation of the electrical conductance of networks. In Sec. 3 a simple phys- 
ical picture gives a theoretical explanation of the results. Section 4 deals with 
the number of link- independent paths as a form of frictionless transport. In 
Section 5 we extend the study of frictionless transport to the case of multi- 
ple sources. Finally, in Sec. 6 we present the conclusions and summarize the 
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results in a coherent picture. 



2 Transport in complex networks 

Most of the work done so far regarding complex networks has concentrated 
on static topological properties or on models for their growth [4,13,11,19]. 
Transport features have not been extensively studied with the exception of 
random walks on specific complex networks [20,21,22]. Transport properties 
are important because they contain information about network function [23]. 
Here we study the electrical conductance G between two nodes A and B of 
Erdos-Renyi and scale-free networks when a potential difference is imposed 
between them. We assume that all the links have equal resistances of unit 
value [24]. 

We construct the Erdos-Renyi and scale-free model networks in a standard 
manner [15,25]. We calculate the conductance G of the network between two 
nodes A and B using the Kirchhoff method [26], where entering and exiting 
potentials are fixed to Va = 1 and Vb = 0. 

First, we analyze the probability density function (pdf) $(G) which comes 
from the probability that two nodes on the network have conductance 

between G and G -\- dG. To this end, we introduce the cumulative distribution 
F{G) = J^^{G')dG', shown in Fig. 2 for the Erdos-Renyi and scale-free 
(A = 2.5 and A = 3.3, with kmin = 2) cases. We use the notation ^sf{G) and 
Fsf(G) for scale-free, and $er(G) and Fer{G) for Erdos-Renyi. The function 
Fsf{G) for both A = 2.5 and 3.3 exhibits a tail region well fit by the power 



5 



law 



(3) 



and the exponent {gc — 1) increases with A. In contrast, Fer(G) decreases 
exponentially with G. 

We next study the origin of the large values of G in scale-free networks and 
obtain an analytical relation between A and go- Larger values of G require 
the presence of many parallel paths, which we hypothesize arise from the high 
degree nodes. Thus, we expect that if either of the degrees kA or ks of the 
entering and exiting nodes is small (e.g. kA > ks), the conductance G between 
A and B is small since there are at most k different parallel branches coming 
out of a node with degree k. Thus, a small value of k imphes a small number of 
possible parallel branches, and therefore a small value of G. To observe large 
G values, it is therefore necessary that both kA and ks he large. 

We test this hypothesis by large scale computer simulations of the conditional 
pdf $sf(G|A;a, ks) for specific values of the entering and exiting node degrees 
kA and ks- Consider first ks <^ kA, and the effect of increasing ks, with kA 
fixed. We find that $sf(G|/ca, /c_b) is narrowly peaked (Fig. 3(a)) so that it is 
well characterized by G*, the value of G when $sf is a maximum. We find 
similar results for Erdos-Renyi networks. Further, for increasing ks, we find 
[Fig. 3(b)] G* increases as G* ~ k%, with a = 0.96 ± 0.05 consistent with the 
possibility that as N ^ oo, a — 1 which we assume henceforth. 

For the case of ks ^ kA, G* increases less fast than ks, as can be seen in Fig. 4 
where we plot G*/kB against the scaled degree x = kA/ks- The collapse of 
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G* /ks for different values of kA and ks indicates that G* scales as 

G- . k.f ["f) . (4) 



Below we study the possible origin of this function. 



3 Transport backbone picture 



The behavior of the scaling function f{x) can be interpreted using the fol- 
lowing simplified "transport backbone" picture [Fig. 4 inset], for which the 
effective conductance G between nodes A and B satisfies 

-- — ^ — ^— , (5) 



G Ga Gtb G 



B 



where 1/Gth is the resistance of the "transport backbone" while 1/Ga (and 
1/Gb) are the resistances of the set of finks near node A (and node B) not 
belonging to the "transport backbone". It is plausible that Ga is linear in kA, 
so we can write Ga — ckA- Since node B is equivalent to node A, we expect 
Gb — cks- Hence 

C ^ I = k (^kA/kB . . 

l/ckA + l/ckB + 1 1 Gtb 1 + kAlkB + ckAlGtb ' ^ ^ 



so the scafing function defined in Eq. (4) is 

•^^^^ " l^X^ckAlGtb ^ T+^' 



The second equafity follows if there are many parallel paths on the "transport 
backbone" so that <C 1/c/c^. The prediction (7) is plotted in Fig. 4 for 
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both scale-free and Erdos-Renyi networks and the agreement with the simu- 
lations supports the approximate validity of the transport backbone picture 
of conductance in scale-free and Erdos-Renyi networks. 

Within this "transport backbone" picture, we can analytically calculate Fsf{G) 
The key insight necessary for this calculation is that G* ~ ks, when ks < 
kA, and we assume that G ~ ks is also valid given the narrow shape of 
$sf(G|A;a, ks)- This implies that the probability of observing conductance G 
is related to ks through ^sF{G)dG ~ M{kB)dkB, where M{kB) is the prob- 
ability that, when nodes A and B are chosen at random, ks is the minimum 
degree. This can be calculated analytically through 

M(A;b)~P(M / P{.kA)dkA (8) 

kB 



Performing the integration we obtain for G < Gmax 

^SF{G)r^G-'^ [gG^2X-l]. (9) 



Hence, for Fsf{G), we have Fsf{G) ~ G~^'^'^~'^\ To test this prediction, wc 
perform simulations for scale-free networks and calculate the values of go — 1 
from the slope of a log- log plot of the cumulative distribution Fsf{G). Prom 
Pig. 5(b) we find that 

c/G - 1 = (1-97 ± 0.04)A - (2.01 ± 0.13). (10) 



Thus, the measured slopes are consistent with the theoretical values predicted 
by Eq. (9). 



4 Number of link-independent paths: transport without friction 

In many systems, it is the nature of the transport process that the particles 
flowing through the network links experience no friction. For example, this is 
the case in an electrical system made of super-conductors [27], or in the case 
of water flow along pipes, if frictional effects are minor. Other examples are 
flow of cars along traffic routes, and perhaps most important, the transport 
of information in communication networks. Common to all these processes 
is that, the quality of the transport is determined by the number of link- 
independent paths leading from the source to the destination (and the capacity 
of each path), and not by the length of each path (as is the case for simple 
electrical conductance). In this section, we focus on non- weighted networks, 
and define the conductance, as the number of link-independent paths between 
a given source and destination (sink) A and B. We name this transport process 
as the max- flow model, and denote the conductance as Gmf- Fast algorithms 
for solving the max-flow problem, given a network and a pair (^4, B) are well 
known within the computer science community [28] . We apply those methods 
to random scale-free and Erdos-Renyi networks, and observe similarities and 
differences from the electrical conductance transport model. Max-flow analysis 
has been applied recently for complex networks in general [18,29], and for 
the Internet in particular [30], where it was used as a signiflcant tool in the 
structural analysis of the underlying network. 

We flnd, that in the max-flow model, just as in the electrical conductance 
case, scale-free networks exhibit a power-law decay of the distribution of con- 
ductances with the same exponent (and thus very high conductance values 
are possible), while in Erdos-Renyi networks, the conductance decays expo- 
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nentially (Fig. 6(a)). In order to better understand this behavior, we plot the 
scaled-flow GuF/ks as a function of the scaled-degree x = kA/ks (Fig. 6(b)). 
It can be seen that the transition at a; = 1 is sharp. For all x < 1 (/c^ < A;^), 
Gmf — X (or Gmf = ^a), while for x > 1 {ks < /ca), Gmf — 1 (or Gmf = ks)- 
In other words, the conductance simply equals the minimum of the degrees of 
A and B. In the symbols of Eq. (6), this also implies that c ^ 1; i.e. scale-free 
networks are optimal for transport in the max-flow sense. The derivation lead- 
ing to Eq. (9) becomes then exact, so that the distribution of conductances is 
given again by $mf,sf(Gmf) ~ G^p^~^\ 

We have so far observed that the max-flow model is quite similar to electrical 
conductance, both have similar flnite probability of flnding very high values of 
conductance. Also, the fact that the minimum degree plays a dominant role 
in the number of link-independent paths makes the scaling behavior of the 
electrical and frictionless problems similar. Only when the conductances are 
studied as a function of distance, some differences between the electrical and 
frictionless cases begin to emerge. In Fig, 7(a), we plot the dependence of the 
average conductance Gmf with respect to the minimum degree min(/c^, ks) of 
the source and sink, for different values of the shortest distance Iab between A 
and B, and flnd that Gmf is independent of £ab as the curves for different Iab 
overlap. This result is a consequence of the frictionless character of the max- 
flow problem. However, when we consider the electrical case, this independence 
disappears. This is illustrated in Fig, 7(b), where G is also plotted against 
the minimum degree mm{kA, ks), but in this case, curves with different iAB 
no longer overlap. From the plot we flnd that G decreases as the distance 
increases. 

To test the validity of our model results in real networks, we measured the 
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conductance Gj^p on the most up to date map of the Autonomous Systems 
(AS) level of the Internet structure [31]. From Fig. 8 we find that the slope 
of the plot, which corresponds to go — 1 from Eq. (9), is approximately 2.3, 
implying that ARi2.15±0.05,in agreement with the value of the degree distri- 
bution exponent for the Internet observed in [31]. Thus, transport properties 
can yield information on the topology of the network. 

5 Multiple Sources 

In many cases it is desirable to explore the more general situation, where the 
transport takes place between groups of nodes, not necessarily between a single 
source and a single sink. This might be a frequent scenario in systems such 
as the Internet, where data should be sent between a group of computers to 
another group, or in transportation network, where we can consider the quality 
of transport between e.g. countries, where the network of direct connections 
(flights, roads, etc.) between cities is already estabhshed. 

The generalization of the above defined models is straightforward. For the 
electrical transport case, once we choose our n sources and n sinks, we simply 
wire the n sources together to the positive potential V — 1, and the n sinks 
to V — 0. For the max flow case, we connect the n sources to a super-source 
with infinite capacity links, the same for the sinks. We then consider the max 
flow between the super-source and the super-sink. 

For simplicity of notation, in this section we denote the electrical conductance 
with G and the max-flow with simply / (instead of Gmf)- 

The interesting quantity to consider in this case is the average conductance. 
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or flow, per source, i.e. G/n or f/n (the averaging notation will be omitted 
when it is clear from the context), since this takes into account the obvious 
increase in transport due to the multiple sources. Considering the transport 
per source gives thus a proper indication of the network utilization. 

In Fig. 9 we analyze the dependence of the transport on n. Consider first the 
case of max flow, in 9(a). It is seen, that the flow first increases with n, and 
then decreases, for both network models and network sizes. This observation 
is consistent with the transport backbone picture described above, if we gen- 
eralize the definition of the transport backbone to include now all nodes in 
the network that are nor sources neither sinks. For small n (or more precisely, 
n <^ N), the transport backbone remains as large, and hence as a good con- 
ductor, as it is in the case of a single source (n = 1). Therefore, we expect 
again the flow to depend only on the degrees of the sources and the sinks. With 
one source and sink, the flow is conflned by the smaller degree of the source 
and the sink. With multiple sources and sinks, the flow will be the minimum of 
[sum of the degrees of the sources, sum of degrees of the sinks] . The more sources 
we have (larger n), the more chances for the sum of the degrees to be larger 
(since for a sum of many degrees to be smaller, we need all the degrees to 
be small, which is less probable). Thus for small n, f/n increases. An exact 
theory for this region is developed below. As n grows larger, not only that 
the transport backbone becomes smaller, but the number of paths that need 
to go through it grows too. Therefore, many paths that were parallel before 
now require the simultaneous usage of the same link, and the backbone is 
no longer a perfect conductor. In other words, the interaction between the 
outgoing paths causes the decrease in the backbone transport capability. 

From these two contradicting trends emerges the appearance of an optimal n, 
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n* . which is the n for which the flow per source is maximized. This has an 
important consequences for networks design, since it tells us that a network 
has an optimal "number of users" for which the utilization of the network is 
maximized. Increasing or decreasing the number of nodes in each group, will 

force each node, on average, to use less of the network resources. 

In Fig. 9(b), we plot the same for the electrical conductance case. But, in 
contrast to the flow case, here the conductance per source only increases with 
respect to the number of sources n. The reason is, that for large n, the number 
of parallel paths decrease, but their lengths decrease (since for so many sources 
and sinks, there is more probability for a direct or almost direct connection 
between some source and some sink), and therefore the conductance of each 
path signiflcantly increase, such that the total flow increases too. All this did 
not affect the flow, since there, as mentioned in Section 4, the total flow do 
not depend on the distance between the source and the sink. 

This fact, that electrical conductance improves with more sources, while flow 
only degrades, actually points out a fundamental difference between the two 
types of transport phenomena. This actually is a consequence of the arguments 
given in Section 4 as for the distance dependence, and completes the picture 
considering the comparison between the models. 

For the max- flow model, a closed form formulas for f{n) and the probability 
distribution of the flow $mf(/) (for a given n) can be derived analytically in 
the region where n <^ N, where we can assume there is no interaction at all 
between the paths of the n pairs, such that the flow is just the minimum of 
[sum of the degrees of the source, sum of the degrees of the sinks] . A compar- 
ison between the theoretical formula and the simulation results will enable us 
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to directly test this hypothesis; moreover it will mark the region {n values) 
where "no interaction" between flow paths exist. 



Next, we present the analytic derivation for small n. Denote by kn the sum 
of degrees of n nodes. The total flow between n sources and n sinks is then 
given by / = min[A;„ (sources), A;„ (sinks)]. To calculate the distribution of flows 
(from which the average flow is found easily) one has first to calculate the 
distribution of the sum of n iidrv's ki. Pki{k) ~ k~'^, [/Cmin ^ k < /Cmax], for SF 
graphs, or Pk^{k) = ^^-'-^^ — for ER graphs. For n » 1, this is easy with the 
central limit theorem. However, for very large n, the effect of the interactions 
will begin to appear. Therefore, one should calculate the distribution of the 
sum directly, using convolution. If we consider the SF degree distributions to 
be discrete. The pdf of a sum of degrees of 2 nodes is then given by - 

h — h - 

PkM= E PkAj) ■ PkAk - j) ; 2kmin<k<2kmax: (H) 

which can be computed exactly in a computer, for a given A^, kmin, A. The cal- 
culation can be easily extended to n equals other pairs of 2. For ER networks, 
it is known that the sum of n Poisson variables with parameter A is also a 
Poisson variable with parameter nX. The probability distribution of the fiow 
itself is readily calculated as the pdf of the minimum of the k^s, analogous to 
Eq. (8)- 



$MF,n(/) = 2-P,„(/) 



J=f 



[PkM)f- (12) 



This can also be computed exactly in a computer, together with the average 
flow (/(n) = I]/=i / ■ ^MF,n(/))- However, for ER networks the sum can be 
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solved: 



^MF,n(/) = 2 ^1 - -YUy- - 2 J\ ) ' ^^^^ 

where r(a) is the Gamma Function and r{a,x) is the Incomplete Gamma 
Function. In Fig. 10, the "no-interactions" theory is compared to the simu- 
lation results, for small n. In panel (a), the pdf of the flows is shown for an 
ER network. The agreement between the theory and simulation is evident. 
In panel (b), we show the average flow per source, f/n, as computed using 
Eq. (12), together with the simulation results, for ER and SF networks. The 
range of applicability of the no-interactions assumptions is clearly seen. 

Scaling law for the conductance with multiple sources, is seen when properly 
scahng the conductance and the number of sources. Plotting G/N vs. n/N, 
(Fig. 11), all curves with different values of N collapse. Thus we conclude the 
scaling form: 

Gr^Ng(n/N), (14) 

where g{x) is a function of a single variable only. The same scaling appears 
for the flow as well (data not shown) . 

6 Summary 

In summary, we find that the conductance of scale-free networks is highly 
heterogeneous, and depends strongly on the degree of the two nodes A and B. 
We also find a power- law tail for ^sf{G) and relate the tail exponent go to the 
exponent A of the degree distribution P{k). This power law behavior makes 
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scale- free networks better for transport. Our work is consistent with a simple 
physical picture of how transport takes place in scale-free and Erdos-Renyi 
networks, which presents the 'transport backbone' as an explanation to the 
fact that the transport virtually depends only on the smallest degree of the 
two nodes. This scenario appears to be valid also for the frictionless transport 
model, as clearly indicated by the similarity in the results. We analyze this 
model and compare its properties to the electrical conductance case. We also 
compare our model results on a real network of the AS Internet and obtain 
good agreement. We then extend the study to transport with multiple sources, 
where the transport takes place between two groups of nodes. We find that this 
mode arises a fundamental difference between the behavior of the electrical 
and frictionless transport models. We also find an optimal number of sources 
in which the transport is most efficient, and explain its origin. 



Finally, we point out that our study can be further extended. We could find 
a closed analytical formula, Eqs. (12), (13) for the distribution of flows in ER 
networks only. It would be useful to find an analytical formula also for SF 
networks, and in addition, to find a analytical form for the average fiow (and 
conductance), per source. In other words, we need to find the scaling function 
g{x) in Eq. (14). Also, it is of interest to investigate the dependence of the 
optimal number of sources, n*, in the network size and connectivity. Finally, 
other, maybe more realistic, models for transport in a network with multiple 
sources should be considered. 
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Fig. 1. (a) Schematic of an Erdos-Renyi network of = 12 and p = 1/6. Note 
that in this example ten nodes have k = 2 connections, and two nodes have k = 1 
connections. This illustrates the fact that for Erdos-Renyi networks, the range of 
values of degree is very narrow, typically close to k. (b) Schematic of a scale-free 
network of N = 12, femin = 2 and A « 2. We note the presence of a hub with 
kmax = 8 which is connected to many of the other nodes of the network. 
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Fig. 2. Comparison for networks with TV = 8000 nodes between the cumulative 
distribution functions of conductance for the Erdos-Renyi and the scale-free cases 
(with A = 2.5 and 3.3). Each curve represents the cumulative distribution F(G) vs. 
G. The simulations have at least 10^ realizations. 
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Conductance G degree of node B 



Fig. 3. (a) The pdf ^sF{G\kA,kB) vs. G for iV = 8000, A = 2.5 and kA = 750 
{kA is close to the typical maximum degree kmax = 800 for N = 8000). (b) Most 
probable values G* , estimated from the maximum of the distributions in (a) , as a 
function of the degree kB- The data support a power law behavior G* ~ feg with 
a = 0.96 ±0.05. 
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Fig. 4. Scaled most probable conductance G* /ks vs. scaled degree x = kA/ks 
for system size AT = 8000 and A = 2.5, for several values of kA and ks' □ 
{kA = 8, 8 < kB < 750), {kA = 16, 16 < fcs < 750), A {kA = 750, 
i < kB < 128), O (fes = 4, 4 < fe^ < 750), y (^s = 256, 256 < kA < 750), 
and > {kB = 500, 4 < fc^ < 128). The curve crossing the symbols is the predicted 
function G* /kB = f{x) = cx/{l + x) obtained from Eq. (7). We also show G* /kB 
vs. scaled degree x = kA/ks for Erdos-Renyi networks with k = 2.92, 4 < fc^ < 11 
and kB = 4: (symbol •). The curve crossing the symbols represents the theoretical 
result according to Eq. (7), and an extension of this line to represent the limiting 
value of G* /kB (dotted-dashed line). The probability of observing fc^ > 11 is ex- 
tremely small in Erdos-Renyi networks, and thus we are unable to obtain significant 
statistics. The scaling function /(x), as seen here, exhibits a crossover from a lin- 
ear behavior to the constant c (c = 0.87 ± 0.02 for scale-free networks, horizontal 
dashed line, and c = 0.55 it 0.01 for Erdos-Renyi, dotted line). The inset shows 
a schematic of the "transport backbone" picture, where the circles labeled A and 
B denote nodes A and B and their associated links which do not belong to the 
"transport backbone" . 
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Fig. 5. (a) Simulation results for the cumulative distribution Fsf(G) for A be- 
tween 2.5 and 3.5, consistent with the power law Fgp Q-(9g-^) (^f.f_ Eq. (9)), 
showing the progressive change of the slope go — (b) The exponent go — I 
from simulations (circles) with 2.5 < A < 4.5; shown also is a least square fit 
qg — ^ = (1.97 lb 0.04) A — (2.01 it 0.13), consistent with the predicted expression 
5G-l = 2A-2 [cf. Eq. (9)]. 
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Fig. 6. (a) Cumulative distribution of link-independent paths (conductance) 
-^mf(Gmf) vs. Gmf compared with the electrical conductance distributions taken 
from Fig. 2. We see that the scaling is indeed the same for both models, but the 
proportionality constant of -Fmf(Gmf) vs. Gmf is larger for the frictionless problem, 
(b) Scaled number of independent paths Gmf/^b as a function of the scaled degree 
kA/kB for scale- free networks of = 8000, A = 2.5 and fcmin = 2. The behavior is 
sharp, and shows how Gmf is a function of only the minimum k. 
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Fig. 7. (a) Average conductance Gmf vs. minimum degree mm(kA, ks) of the source 
and sink A and B for different values of the shortest distance Iab- The relation 
is independent of ^ab indicating the independence of Gmf on the distance. The 
network has = 8000, A = 2.5, fcmin = 2. (b) Average conductance G vs. minimum 
degree mm{kA, ks) of the source and sink A and B for different values of distance 
Iab- The independence of G with respect to Iab breaks down and, as Iab increases, 
G decreases. Once again, N = 8000 and A = 2.5, but the average has been performed 
for various ks < kA and kA = 750. 
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Fig. 8. Cumulative distribution F(G^p) vs G^^p for the Internet. This slope is in 
good agreement with the scale- free structure that has been observed for the Internet 
(see text). 
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Number of sources n number of sources n 

Fig. 9. (a) The average flow per source, f jn vs. the number of sources n, for ER and 
SF networks, and two values of network size N . For all curves there is an optimal 
n*, for which the flow per source is maximized. Obviously, n* grows with such 
that the larger the network is, more users can use it in an optimal way. Above n*, 
the flow per source decreases, (b) The average conductance per source, Gjn vs. the 
number of source n again for ER, SF networks, and different network sizes. Here 
there is no optimal point- The more users in the network the more useful it is for 
transport. See the text for the qualitative explanation of this behavior. 
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Fig. 10. (a) For an ER network, with average degree k = 4, N = 8192, we plot 
the probability distribution of the flows <l*MF,m for given n = 4,8, 16, as wefl as the 
theoretical prediction from Eq. (13). (b) For ER network with k = 8, N = 4096, and 
SF network with A = 2.4, k^i^ = 2, N = 8192, we plot the average flow per source, 
/ vs. n, and compare it with the theoretical prediction under the assumptions 
of no interactions between the parallel paths. The number of sources where this 
assumption no longer holds is marked. 
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Fig. 11. (a) The average conductance G is plotted vs. the number of sources n for 
ER networks with average degree k=A, and various values of network sizes N. (b) 
The same data, plotted now as G/N with respect to n/N. With this scaling all the 
curves collapse, implying the scaling law G ^ Ng{n/N). 
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